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With the advent of imaging atmospheric Cherenkov telescopes in late 1980's, ground-based 
observation of TeV gamma-rays came into reality after struggling trials by pioneers for twenty 
years, and the number of gamma-ray sources detected at TeV energies has increased to be over 
seventy now. In this review, recent findings from ground-based very-high-energy gamma-ray 
observations are summarized (as of 2008 March), and up-to-date problems in this research field 
are presented. 
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1. Introduction 

There has been a dramatic progress in very-high- 
energy-*^^ gamma-ray astronomy in the last decade. With 
the advent of imaging atmospheric Cherenkov telescopes, 
the number of gamma-ray sources detected in the TeV 
energy range has increased significantly in the last sev- 
eral years, and now it exceeds seventy, according to the 
rapporteur talk by Jim Hinton at the International Cos- 
mic Ray Conference in 2007.^^ These gamma-ray sources 
form a new class of high-energy objects in the Universe, 
including active galactic nuclei, radio galaxies, galactic 
binaries, pulsar wind nebulae, in addition to supernova 
remnants which are assumed to be the origin of cosmic- 
rays for a long time. Exploring the emission mechanism 
from these objects is a big challenge in astrophysics. 
Non-thermal nature of emission inherently needs multi- 
wavelength observations to study the phenomena, involv- 
ing astronomers working in other wavelength. Further- 
more, many unidentified gamma-ray sources have been 
reported and are posing new mysteries. 

In this review, highlights of recent findings from 
ground-based very-high-energy gamma-ray observations 
are summarized, and up-to-date problems in this re- 
search field are presented. We concentrate on recent re- 
sults from Cherenkov telescopes here, but new results 
from particle shower arrays are also exciting (G. Sinnis 
and M. Takita, in these proceedings.) Recent review on 
this field is also found elsewhere.'^-' 

2. Imaging Atmospheric Cherenkov Telescopes 

Ground-based imaging Cherenkov telescopes are be- 
coming a powerful tool to study very high energy gamma- 
rays with their capability to discriminate gamma-rays 
from background hadrons (protons and nuclei).^'' 

Gamma-ray images come from purely electromagnetic 
showers and are sharp and oriented toward the ob- 
ject being tracked. They can be separated from cosmic- 
ray showers using imaging parameters first defined by 
A.M. Hillas:^^ width, length, distance, asymmetry. Distri- 
butions of theses parameters for gamma-ray showers are 
different from those for hadronic showers, and we can ex- 
tract gamma-ray signals statistically from observed im- 
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ages. The first firm detection of a TeV gamma-ray signal 
from the Crab nebula by the Whipple group^-' utilized 
this imaging technique. 

Stereoscopic observation of Cherenkov images came 
into practical use by the HEGRA group. Incoming di- 
rection of a gamma-ray can be determined by intersec- 
tion of axes of elongated Cherenkov images observed by 
multiple telescopes, separated by about 100 m, which is 
the size of light pool of Cherenkov light fiash, more pre- 
cisely compared to single telescope observation. Distance 
of centroid of shower images from the incoming direction 
is a measure of shower maximum height in the atmo- 
sphere, which also allows better estimation of gamma- 
ray energy. In addition, background Cherenkov hashes 
caused by local cosmic-ray muons traversing near tele- 
scopes can be effectively rejected requiring coincidence 
between telescopes. 

Table I summarizes characteristics of Cherenkov tele- 
scope systems in operation. 

3. Recent Topics 

Table II is a summary of TeV sources which are fairly 
established (detected by multiple groups and/or detected 
at high significance). Over 40 sources in this table were 
found by the H.E.S.S. group, mostly by their Galactic 
plane survey.®^ Note the classification of sources are not 
unique and depend on person who makes a list: identifi- 
cation of TeV sources with objects found in other wave- 
length are sometimes difficult, because of the limited 
angular resolution of TeV observation: therefore many 
TeV sources are left unidentified. (But there are 'dark' 
sources: see section 3.1.6.) 

Fig. 1 is a skymap of TeV sources in the Galactic co- 
ordinates compiled by R. Wagner.^"' Concentration of 
sources along the Galactic plane is obvious, but please 
note the exposure time is far from uniform. 

3.1 Galactic Sources 
3.1.1 Supernova Remnants 

Needless to say, supernova remnants (SNRs) are long 
considered to be primary sources for galactic cosmic 
rays.^^-* They are energetic enough to support the cos- 
mic ray luminosity of the Galaxy and their sizes are 
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Table I. Current Cherenkov telescope systems 



Group 


Location 


Latitude 


Longitude 


Height 


Telescopes 


Start 


SHALON 


Russia 


42° S 


75°E 


3,338 m 


llm^x2 


1992 


TACTIC 


India 


25° S 


73° E 


1,300 m 


9.5m^x4 


2000 


CANGAROO-III 


Australia 


31°S 


137° E 


160 m 


57m^x4 


2004 


H.E.S.S. 


Namibia 


23° S 


16.5°E 


1,800 m 


107m^x4 


2004 


MAGIC 


Canary Is. 


29°N 


18°W 


2,200 in 


237 X 1 


2004 


VERITAS 


Arizona 


32° N 


111°W 


1,300 m 


110in^x4 


2007 




Fig. 1. TeV skymap (circa early 2008) compiled by R. Wagner.^''' 



Table II. Summary of TeV sources as of March 2008^^^ 



Source type 


Count 


Galactic sources 




Shell-type supernova remnant 


3 


Supernova remnant 


11 


Pulsar wind nebula 


16 


Binary 


4 


Galactic center 


1 


Galactic ridge 


1 


Young stellar cluster 


1 


Unidentified 


14 


Extragalactic sources 




High-frequency peaked BL Lac (HBL) 


17 


Low-frequency peaked BL Lac (LBL) 


1 


Radio galaxy 


1 


Flat-spectrum radio quasar (FSRQ) 


1 


Total 


71 



large enough to confine cosmic rays during acceleration 
by diffusive Shockwave. Bright synchrotron radio emis- 
sion from many SNRs indicate existence of non-thermal, 
high-energy electrons. Evidences of acceleration of pro- 
tons in SNR, however, were rare. 

The first SNR detected at TeV energies is RX 
J1713.7-3946 (by CANGAR00-I,i2) CANGAROO- 
11,^^) and H.E.S.S.^^)) where the shell structure was 
clearly observed. Although the spectrum is better 
explained by emission from neutral pions produced in 
hadronic interaction of protons, ^^"^^^ inverse Compton 
emission by electrons is still a viable process to ac- 



count for the TeV spectrum. RX J0852.0-4622 (by 
CANGAROO-II^s) and H.E.S.S.^^)) and RCW 86 (by 

H.E.S.S."'^*^) arc other SNRs which show shell-like emis- 
sion profile, well correlated with X-ray intensity profiles. 
Cas A, first claimed by HEGRA^^^ and confirmed to 
be a TeV emitter by MAGIC,^") and CTA 37A/37B, 
recently claimed by H.E.S.S. could also be shell-like 
SNRs, but their small angular sizes do now allow fur- 
ther study. These shell-like structure arc often well cor- 
related with X-ray images, and one may suppose TeV 
emission comes from inverse Compton scattering of low 
energy photons by high-energy electrons which radiate 
synchrotron X-rays. Still, gamma-rays derived from 
neutral pion decays produced in collision of high-energy 
protons with ambient matter could account for TeV emis- 
sion, since amplified magnetic fields in SNRs might 
suppress inverse Compton emission. 

Uchiyama et al. analyzed the X-ray images of RX 
J1713. 7-3946 observed by Chandra in multiple years, and 
found there are 'hot spots' showing time variation in one- 
year scale. ^■''^ This time scale, compared with the cooling 
time of high energy electrons, means the ambient mag- 
netic field is stronger than 1 mG. This magnetic field sup- 
presses the inverse Compton emission in the TeV region, 
and may indicate TeV emission is produced by hadronic 
process. However, Butt et al. argue that the existing mul- 
tiwavelength data on this object do not support such 
conclusion. 

On the other hand, TeV emission profiles of IC443 (by 
MAGIC27) and VERITAS28)) andW28 (by H.E.S.S.^^)) 
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Tabic III. TcV sources associated with pulsar wind nebulae.^' ®^ 
These associations were established through a range of methods, 
which are given in the table in abbreviated form: Pos: The posi- 
tion of the centroid of the VHE emission can be established with 
sufficient precision that there is no ambiguity as to the low en- 
ergy counterpart. Mor. There is a match between the gamma-ray 
morphology and that seen at other (usually X-ray) wavelengths. 
EDMor. Energy-dependent morphology which approaches the 
position/morphology seen at other wavelengths at some limit, 
and is consistent with our physical understanding of the source. 



Object 


Pulsar 


Method 


Discovered 


Crab nebula 


B0531-I-21 


Pos 


Whipple, 1989 


MSH 15-52 


B1509-58 


Mor 


HESS, 2005 


VelaX 


B0833-45 


Mor 


HESS, 2006 


HESS J1825-137 


B1823-13 


EDMor 


HESS, 2005 


PSR J1420-6049 


J1420-6048 


Mor 


HESS, 2006 


The Rabbit 


J1420-6048 


Mor 


HESS, 2006 


G0.9-f0.1 




Pos 


HESS, 2005 



do not show shell-like structure and seems to correlate 
with molecular cloud profiles observed by a CO emission 
line. This interaction may indicate the hadronic origin 
of gamma-rays producicd in collisions of high-energy pro- 
tons with molecular cloud as targets. 

HESS J0632+057 was found close to the rim of 
the Monoccros Loop SNR/Rosette Nebula region. '^"^ 
It is point-like and has no clear counterpart at other 
wavelengths, but is possibly associated with the GeV 
source 3EG J0634-I-0521, a weak X-ray source IRXS 
J063258.3+054857 and the Be-star MWC 148. 

Thus, although supernova remnants are proved to be 
the sites producing high energy particles in the Universe, 
species of accelerated particles are still unclear:^^) solving 
the long-standing mystery of cosmic ray origin is still an 
important issue of gamma-ray astronomy. 

3.1.2 Pulsar Wind Nebulae 

Table III is a list of TeV sources with well established 
association with pulsar wind nebulae. Including associa- 
tions with weaker evidence, 18 out of 71 sources detected 
at TeV energies seem to be associated with pulsar wind 
nebulae, which was a rather surprising discovery re- 
vealed by the H.E.S.S. Galactic survey.^) Most of them 
are associated with relatively young (< 10^ years) and 
large spin-down pulsars, "^^^ which means the gamma-ray 
luminosity is supported by the pulsar spin-down energy. 
Their profiles show extended (order of 10 pc) structure, 
often displaced from pulsar positions. No pulsation has 
been reported in the TeV region, even for the Crab neb- 
ula which is a strong TeV source and best studied by 
a number of groups. Thus the TeV emission are nat- 
urally ascribed to inverse Compton emission by high- 
energy electrons accelerated in the vicinity of pulsars. 

This hypothesis on emission mechanism is further sup- 
ported by the energy-dependent morphology in HESS 
J1825-137.^^^ The photon indices from a power-law fit 
in different regions show a softening of the spectrum with 
increasing distance from the pulsar. The observed energy 
dependent morphology may be an evidence for cooling of 
electrons in the nebula. 

HESS J1837— 069, which was classified as unidenti- 
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fied before, has now been added to this category with 
the discovery of an energetic 70.5 ms pulsar in AX 
J1838.0-0655 using X-ray data by RXTE.^^) 

3.1.3 Gamm,a-ray Binaries 

PSR B1259-63/SS2883 This is a 48 ms pulsar and 
a Be star binary in a highly eccentric orbit. Gamma-ray 
flux at detectable level was predicted when the binary is 
near the periastron passage via interaction of pulsar wind 
with the radiative environment of the binary system.^^^ 
H.E.S.S. detected such a variable flux around the 2004 
periastron,^^) with a double-peak light curve as predicted 
by some models. ^''^ 

LS5039 This is a high mass X-ray binary compris- 
ing a massive star and a compact object, and is resolved 
into a bipolar radio outflow emanating from a central 
core, thus it is often classified as the microquasar class. 
H.E.S.S. detected a modulated gamma-ray signal with 
the 3.9 day orbital period in 2005. The emission max- 
imum {(/) ~ 0.7) appears to lag behind the apastron 
epoch and to align better with the inferior conjunction 
{(f) = 0.716), when the compact object lies in front of 
the massive star. The flux minimum occurs at phase 
{(j) ~ 0.2), slightly further along the orbit than supe- 
rior conjunction (0 = 0.058). This behavior is not easy 
to explain by a simple assumption and offers a challenge 
to model builders. 

LSI+61 303 MAGIC observed this microquasar, 
known in the GeV region as 2CG 135+01 (by COS B) 
or 3EG J0241-h6103 (by EGRET), for six orbital cycles 
(forb = 26.5 days) and obtained a modulated gamma-ray 
signal which is significant only between orbital phase 0.4 
and 0.7.39) j.^^ 

maximum is detected at phases 0.5- 
0.6, overlapping with an X-ray outburst and the onset 
of a radio outburst, is shifted from the phase when the 
two stars are closest to one another, implying a strong 
orbital modulation of the emission or the absorption pro- 
cesses. The VERITAS group confirmed the modulated 
emission. 

Cyg X-1 This well-known black hole X-ray binary 
was observed by MAGIC for a total of 40 hours dur- 
ing 26 nights in 2006. Although the steady TeV emission 
was not observed, a 4cr-level (3.2cr after trial correction) 
evidence was obtained for one night on September 24, 
coinciding with an X-ray flare. ''^^ The significance of de- 
tection may not be strong enough to claim Cyg X-1 is a 
TeV gamma-ray source.^) 

3.1.4 Stellar Cluster 

The H.E.S.S. group reported Westerlund 2. a young 
open stellar cluster which contains dozen 0-stars and 
two Wolf-Rayet stars, is an extended gamma-ray source, 
HESS J1023-575.''2) This is a new source class and the 
gamma-ray emission mechanism is a new challenge to 
theorists: extended profile suggests there might be col- 
lective effects, possibly of stellar winds in the cluster. 

3.1.5 Galactic Center 

The Galactic Center is a confirmed TeV source (by 
CANGAR00-II,''3) Whipple,'*^) and H.E.S.S.^^)), but 
the nature of the emission is not known. The energy spec- 
trum shows a simple power-law shape, not compatible 
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with dark matter annihilation signal,'*^-' but the position 
coincides with the central radio source, Sgr A*.^'') 

Moreover, when the Galactic center source and the 
nearby source, GO.9+0.1, are subtracted from the 
H.E.S.S. data, emission extended in galactic longitude 
for roughly 2° and also in galactic latitude with a char- 
acteristic width of about 0.2° were secn.'**^ The power- 
law index of the spectrum (2.3) of this Galactic ridge 
gamma-ray emission is harder than the local cosmic-ray 
spectrum (2.7): this may indicate the propagation effects 
are less pronounced than in the Galaxy as a whole due 
to the proximity of particle accelerators. 

3.1.6 Unidentified Sources 

Significant number of the TeV sources distributed near 
the Galactic plane are unidentified. Some have no com- 
pelling counterparts, but some are completely dark in 
other wavelength.^^' Former sources could be identified 
if the angular resolution of Chcrcnkov telescopes arc im- 
proved in future detectors, but latter sources bring up a 
new mystery in astrophysics. Multiwavelength observa- 
tions of these sources are going on to reveal their identi- 
ties. 

HESS J1908+063 was discovered during the ex- 
tended H.E.S.S. survey of the Galactic plane, and it coin- 
cides with the recently reported MILAGRO unidentified 
source MGRO J 1908-1-06. 

TeV J2032+4130 was discovered by HEGRA^i) in 
the Cygnus complex region from their Galactic plane 
survey, and confirmed by Whipple^^^ and MAGIC.^^^ It 
is also within the extended MILAGRO source MGRO 
J2031-t-41.54) 

Funk et al. studied the TeV-GeV connection in Galac- 
tic sources comparing the TeV sources with the EGRET 
catalog. Surprisingly, few common sources are found in 
terms of positional coincidence and spectral consistency. 
Distribution of integrated energy flux (Fig. 2) shows 
almost separate peaks for TeV sources and EGRET 
sources, which means the current TeV sensitivity is much 
better than GcV sensitivity. The TeV upper limits put 
strong constraints on simple power-law extrapolation of 
several of the EGRET spectra and thus strongly suggest 
cutoffs in the unexplored energy range from 10 GeV to 
100 GeV. 

3.2 Extragalactic Sources 

Table IV is a list of known extragalactic sources in the 
TeV region in the order of redshift (z). 

3.2.1 Blazars 

Most of sources listed in Table IV are blazar-type ac- 
tive galactic nuclei (AGN), including the second and 
third established TeV sources, Mrk 421 and Mrk 501. 
The prominent feature of gamma-ray emission from this 
class is the rapid time variability, which should be re- 
lated to the central engine of activity, accretion process 
around supermassive black holes. In many cases, gamma- 
ray emission is accounted by inverse Compton process 
by high-energy electrons, which naturally explains X-ray 
and TeV correlated time variation, but hadronic process 
is still a viable option. 




log,(,(Energy Flux (ergs cm s ^)) 

Fig. 2. Distribution of integrated energy flux for sources in the 
Inner Galaxy. For EGRET the energy flux between 1 GeV and 
10 GeV, for the H.E.S.S. sources, the energy flux between 1 TeV 
and 10 TeV is shown. Also shown is the sensitivity prediction 
for the GLAST-LAT for a typical location in the Inner Galaxy 
{£= 10°,6 = 0°). 



Table IV. Known extragalactic TeV sources as of March 2008.^^ 
(See Table II for abbreviation of source class.) 



Object 


z 


Class 


Discovered 


M87 


0.004 


Radio 


HEGRA, 2003 


Mrk 421 


0.031 


HBL 


Whipple, 1992 


Mrk 501 


0.034 


HBL 


Whipple, 1996 


lES 2344-1-514 


0.044 


HBL 


Whipple, 1998 


Mrk 180 


0.046 


HBL 


MAGIC, 2006 


lES 1959-1-650 


0.047 


HBL 


7TA, 1999 


BL Lac 


0.069 


LBL 


MAGIC, 2007 


PKS 0548-322 


0.069 


HBL 


HESS, 2007 


PKS 2005-489 


0.071 


HBL 


HESS, 2005 


PKS 2155-304 


0.116 


HBL 


Durham, 1999 


H 14264-428 


0.129 


HBL 


Whipple, 2002 


lES 0229-1-428 


0.140 


HBL 


HESS, 2007 


H 2356-309 


0.165 


HBL 


HESS, 2006 


lES 1218-1-304 


0.182 


HBL 


MAGIC, 2006 


lES 1101-232 


0.186 


HBL 


HESS, 2006 


lES 0347-121 


0.188 


HBL 


HESS, 2007 


lES 1011-1-496 


0.212 


HBL 


MAGIC, 2007 


PO 1553+113 


> 0.25 


IIBL 


HESS. 2006 


3(,'279 


0.53(> 




mac;k:, 2007 



In July 2006, the H.E.S.S. group reported an excep- 
tionally strong flare of PKS 2155-304.^^) The peak 
flux was seven times the flux observed from the Crab 
nebula, and the flux varied at time scale as short as one 
minute. The CANGAROO-III group confirmed this 
flaring activity but in a different time zone.^^-* This short 
time variability set limit on the size and Doppler factor 
of gamma-ray emission region near the central massive 
black hole of the blazar. Assuming the emission region 
has a size comparable to the Schwarzschild radius of a 
^ lO'^AfQ black hole, Doppler factors greater than 100 
are required to accommodate the observed variability 
time scales. 

The MAGIC group reported a detection of 3C279,-'^'') 
a prominent flat-spectrum-radio-quasar-type AGN in the 
GeV region observed by EGRET.*^°) It was believed that 
its distance {z = 0.536) was too far to be detected 
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Fig. 3. Observed spectral indices for AGN in the VHE region.^) 

in very-high-energy gamma-rays due to the absorption 
effect by the extragalactic infrared background radia- 
tion (see next subsection) , and we should draw attention 
to further observations of this object, since the current 
statistics is too limited to draw spectral information. 

3.2.2 Extragalactic background radiation 

Stcckcr et al.^^^ pointed out that extragalactic dif- 
fuse infrared background radiation blocks the propaga- 
tion of TeV gamma-ray over large distances (2; ^ 0.1 or 
d ^ 100 Mpc) by producing electron-positron pairs. How- 
ever, measuring the intensity of the background infrared 
radiation is a difficult task, since subtracting foreground 
sources from observed data is a very complicated process. 
Recent progress has been summarized by many authors: 
see ref.,^^^ for example. In turn, we can infer the intensity 
from observation of TeV spectra of distant sources, as- 
suming the spectra at sources are known: this cannot be 
known, in fact, so what we can do in practice is to place 
an upper limit of the intensity of infrared background ra- 
diation. Such studies (ref.,^^' for example) indicate the 
infrared background radiation might be at lower level 
than anticipated. This is a result placing constraint on 
galaxy formation theory, and is a good example of po- 
tential power of TeV gamma-ray astronomy toward cos- 
mology. ^''^ 

Fig. 3 is a plot of spectral indices of AGN in the TeV 
energy region against redshift. If the intrinsic spectra 
of AGN are similar, there should be a rising tendency in 
this plot as redshift increases, but we cannot draw clear 
conclusion from this plot: surely we need more samples. 

3.2.3 M87 

HEGRA reported a detection of M87, a nearby radio 
galaxy in the Virgo cluster (z = 0.00436 or ^ 16 Mpc) 
in 2003.®^-' H.E.S.S. confirmed the gamma-ray emission 
of this source, which showed the short-time variability 
of 2-day scale based on observation between 2003 and 
2006.^^' This time scale is the order of the light crossing 
time of the central black hole. The VERITAS group also 
detected a gamma-ray signal in 2007. 



3.2.4 Centaurus A 

Recent data from the Pierre Auger Observatory sug- 
gests this nearby giant radio galaxy (z — 0.0018) could 
be a source of ultra-high-energy cosmic rays.^^^ If true, 
one may expect it to emit gamma-rays. However, in 
the TeV region, only upper limits on gamma-ray flux are 
reported except one evidence reported in 1970's.™~^^' 

3.2.5 Clusters of Galaxies 

As the largest systems in the Universe, clusters of 
galaxies can harbor high-energy cosmic-rays for cosmo- 
logical time which may be accelerated in merger shocks 
and/or accretion shocks, and can emit gamma-rays at 
detectable level via various possible processes. 

The Whipple group observed Perseus (z = 0.018) and 
Abell 2029 {z = 0.078), but extended gamma-ray emis- 
sion signals were not observed. The H.E.S.S. group 
looked for gamma-ray emission from two nearby clusters, 
Abell 496 {z = 0.033) and Coma (z = 0.024). '^^^ The 
CANGAROO-lII group selected Abell 4038 (z = 0.028) 
and Abell 3667 (z = 0.055)^'") for target. Again, both 
groups could not see a hint of signal of gamma-ray emis- 
sion. 

4. Future Projects 

With these impressive discoveries, next-generation 
large-scale atmospheric Cherenkov telescope systems are 
under hot discussion in recent years to enhance the sen- 
sitivity in TeV gamma-rays furthermore. The H.E.S.S. 
group is constructing a huge 28 m diameter telescope, 
called H.E.S.S. 11, at the center of the present array, 
and the MAGIC II telescope, an advanced copy of the 
present MAGIC 17 m telescope 80 m apart, will be com- 
pleted this year in order to enjoy the merit of stereo- 
scopic observation.''*) In the longer range, the size of 
next-generation projects naturally requires international 
collaboration. CTA (Cherenkov Telescope Array)™) is 
one of such projects and is lead mainly by H.E.S.S. 
and MAGIC collaborators. AGIS (Advanced Gamma- 
ray Imaging System)*^) is another initiative proposed 
mainly by VERITAS collaborators. These projects aim 
to increase the sensitivity by a factor of ten, as well as 
to enhance the gamma-ray energy range both lower and 
higher regions. The trend seen in the Kifune plot (Fig. 4) 
predicts a bright future of TeV astronomy with an order 
of one thousand sources in mid 2010's. 

5. Summary 

Now the very-high-energy window of the Universe is 
open, and Cherenkov telescopes are providing additional 
2-3 decades to the photon spectra of high-energy objects 
in the sky. We can see wider variety of source classes 
among seventy TeV sources than anticipated, indicating 
cosmic accelerators are ubiquitous in the Universe. How- 
ever, understanding the underlying physical processes 
is still at a preliminary stage. Together with GLAST 
(P. Michelson, in these proceedings.), to be launched in 
2008 to observe GeV gamma-rays from space, we can ex- 
pect many findings from the exploration of high-energy 
Universe with next-generation ground-based gamma-ray 
detectors. 
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Fig. 4. Yearly increase of X-ray, gaimna-ray and vcry-high-cncrgy 
(VHE) gamma-ray sources. (This type of plot is called as Kifune 
plot, named after Prof. Tadashi Kifune who invented it.) 
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In this article, very high energy (VHE) means the energies 
from 100 GeV to 10 TeV, and has the same meaning of TeV 
region. 
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